Polyelectrolyte Solutions with Multivalent Salts 
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We investigate the thermodynamic properties of a polyelectrolyte solution in a presence of multi- 
valent salts. The polyions are modeled as rigid cylinders with the charge distributed uniformly along 
the major axis. The solution, besides the polyions, contain monovalent and divalent counterions as 
well as monovalent coions. The strong electrostatic attraction existing between the polyions and 
the counterions results in formation of clusters consisting of one polyion and a number of associated 
monovalent and divalent counterions. The theory presented in the paper allows us to explicitly 
construct the Helmholtz free energy of a polyelectrolyte solution. The characteristic cluster size, as 
well as any other thermodynamic property can then be determined by an appropriate operation on 
the free energy. 
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1. Introduction 

The polyelectrolytes are ubiquitous in our every day 
life. They are present in a variety of industrial prod- 
ucts ranging from the hair conditioners to the superab- 
sorbants used in baby diapers. The carrier of genetic 
code, the DNA, is also a polyelectrolyte. Not withstand- 
ing their omnipresence, our understanding of this class of 
polymers is far from complete. The long range Coulomb 
interaction combined with a strong asymmetry existing 
between the various entities of a polyelectrolyte solution, 
has proven to be a formidable challenge to the standard 
methods of liquid state theory. The scaling approaches 
which were so successful in elucidating the properties of 
non-ionic polymers have, so far, failed in their extension 
to polyelectrolytes. What seems to be lacking is a mean- 
field theory similar to the one constructed by Flory for 
the simple polymers and by Debye and Hiickel for the 
symmetric electrolytes. Recently, we have attempted to 
construct such a mean-field theory for the special class of 
rigid polyelectrolytes, an example of which is provided by 
a solution of DNA segments. The advantage of working 
with a rod-like polyelectrolytes is that the intramolecu- 
lar degrees of freedom associated with the conformational 
state of a polyion can be neglected. The polyelectrolyte 
solution then becomes a strongly asymmetric electrolyte, 
to which a generalization of the Debye-Huckcl-Bjcrrum 
(DHBj) theory can be applied. Additional advantage 
in working with rigid polyelectrolytes is that they obey 
the Manning infinite dilution limiting laws These 
laws are extremely useful since they provide the bound- 



ary conditions that any theory of rod-like polyelectrolytes 
must satisfy. 



We have demonstrated that the extension of the DHBj 
theory to rigid polyelectrolytes, indeed, satisfies the Man- 
ning Limiting Laws Furthermore, the theory is in 
excellent agreement with the experiments on rigid poly- 
electrolyte solutions in a presence of monovalent salt . 
Recently we have extended our work to study the inter- 
action of a polyelectrolyte with a cationic surfactants [EJ . 
The latter problem is of great practical interest in a va- 
riety of biochemical and biomedical applications related 
to design of gene delivery systems ||. Again, the DHBj 
theory has proven successful in giving quantitative agree- 
ment with experiments without a need for any adjustable 
parameters. In the current work we shall extend the the- 
ory to allow for the presence of both the monovalent and 
the divalent salts. This problem, once again, is moti- 
vated by the biomedical applications to gene therapy. It 
has been recently found that introduction of a DNA into 
host cells is facilitated by the presence of a small amount 
of divalent salt such as calcium phosphate || . 



Since the bulk of the calculations is similar to those 
presented in the earlier work, we shall not enter in great 
detail and instead refer the reader to the previously pub- 
lished papers |^-|| . 
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FIG. 1. The DNA binding isotherms for the monovalent 
and divalent counterions. The parameters are: Z = 440, 
i = 4.17, p p = 2 x 10" 6 M, pu = 10" 2 M, a p = 27°A, 
a c = 7.041. 



2. The model 



We begin by defining the Primitive Model of Poly- 
electrolyte (PMP) H . Our system consist of long cylin- 
dric polyions, spherical counterions, monovalent salt, and 
divalent salt, inside the volume V. The polyions, of 
length L, have diameter a p and charge — Zq uniformly 
distributed along the length of the cylinder. The distance 
between the charged groups is b = L/Z. The counterions 
and coions are modeled as rigid spheres of diameter a c 
and charge +q, +2q, —q located at their centers. We shall 
not distinguish between the monovalent counterions de- 
rived from the polymers or salt, and assume that all the 
salt is dissociated. Similarly, to simplify calculations, we 
shall treat the coions derived from the monovalent and 
the divalent salts as being identical. The solvent is mod- 
eled as a continuum medium of dielectric constant D. 
The distance of closest approach, a, between a polyion 
and a counterion is a — (a p + a c )/2. This represents the 
"exclusion cylinder" due to the hard-core repulsion. 

The strong electrostatic attraction existing between 
the macromolecules and the counterions results in some 
number of counterions becoming associated with the 
polyions, producing clusters composed of one polyion, 
niB monovalent counterions, and ds divalent counteri- 
ons. In fact, we expect that there shall be a distribution 
of cluster sizes. Motivated, however, by our previous 
work [Q] we know that the polydispersity in cluster sizes 
is not very important, if one is interested in thermody- 
namics of a polyelectrolyte solution 0. Since not all of 
the counterions condense onto the polyions, for any non- 
zero temperature some free, unassociated, counterions re- 
main in the solution. We denote the density of unassoci- 
ated monovalent counterions as p + = (Z — ms)p P + Pm, 



and the density of unassociated divalent counterions as 
p++ = pd — oIbPp, where p p , pu, and po are the densities 
of the polyions, the monovalent salt, and the divalent salt 
respectively. Since the coions do not participate in asso- 
ciation their density remains unchanged, p_ = pm + ^Pd- 
The goal of the theory is to determine the characteristic 
cluster size, i.e. the values of tub and ds- 

Complete thermodynamic information about the sys- 
tem is contained in its Helmholtz free energy. The con- 
dition that the free energy must be minimum allows 
us to determine the values of tub and ds- Once this 
is done, all the thermodynamic functions of the sys- 
tem can be found through the appropriate operations 
on the free energy. For example, the pressure inside 
the polyelectrolyte solution is a Legendre transform of 
the Helmholtz free energy density / = — F/V, given by 
p(T, {p t }) = f(T, {pt}) + J2 PtPt, where the chemical po- 
tential of a specie of type t (clusters, counterions, and 
coions) is p t = — df/dpt- 

The free energy cannot be calculated exactly. We 
shall, therefore, attempt to construct it as a sum of the 
most relevant contributions. These can be divided into 
an electrostatic and an entropic ones. The electrostatic 
contribution arises as the result of polyion-counterion, 
polyion-polyion, and the counterion-coion interactions. 
The entropic contribution is the result of mixing of vari- 
ous species [[|. 

The free energy for the polyion-counterion interaction 
is calculated in the framework of the Debye-Hiickel (DH) 
theory ||. The full calculation is presented in Ref. 
and here we shall just quote the result, 



(3r = p pi Z- mB -2dBf^{^: 
x < 2 In [KaKi(Ka)] — I(ko) + 
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(1) 



with na = ^f(Airp\/T*), T* = Dk B Ta/q 2 is the reduced 
temperature, pi = pia 3 , and 



I(ko) = 



dx 



K*{x) 



(2) 



where K n is the Bessel function of order n. The modi- 
fication of replacing the total microion density by twice 
the ionic strength, p\ = p + + Ap ++ + p_ , is necessary 
when discussing multivalent salts. 

The contribution to the free energy due to the inter- 
actions among the various free ions is given by the usual 
DH theory g, 



47T ai \ 2 



(3) 



The free, unassociated, counterions and coions screen 
the polyion-polyion interactions, producing an effective 
potential of short range. The polyion-polyion contribu- 
tion to the free energy can then be calculated using a 
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Van der Waals type of approach Q. To this end the 
polyion-polyion contribution is expressed as a second 
virial term, averaged over the relative angle sustained 
by two molecules. The result is, 

= -27ra 3 exp(-2/ta)(Z - m B - 2d B ) 2 p p 2 

The entropic contribution is found using the Flory the- 
ory @- Thus, the increase in free energy due to mixing 
of various species is the sum of their ideal free energies, 
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salt, the condensation is dominated by the monovalent 
ions, since the divalent ions gain more entropy, and thus 
lower the total free energy, by staying free. However as 
the concentration of divalent salt increases, the gain in 
electrostatic energy due to condensation wins over the 
entropy and more divalent ions become associated with 
the macromolecules. We note that unlike the associa- 
tion with the ionic surfactants, which exhibits a large 
degree of cooperativity characterized by the sharp rise in 
the surfactant binding fraction, the replacement of the 
condensed monovalent counterions by the divalent coun- 
terions proceeds quite smoothly. 



where, in obvious notation, the respective volume frac- 
tions are, 
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Q is the internal partition function of an isolated specie 
t. For structureless particles, (+ = C- = C++ = 1- The 
calculation of £ c /, proceeds along the lines given in Ref. 

We find that to a good accuracy the internal par- 
tition function of a cluster can be approximated by 

In £cZ = — £ S (— 2?7ii — 4m2 + 4toiTO2 + m\ + 4:m 2 ) + 
— Z [(1 — mi — 7712) In (1 — mi — 7712) + 
+mi In mi + m.2 In 7712] , (10) 

where the binding fractions are, mi = ttib/Z, iri2 — 
(1b jZ, and S is 



S 
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Z[i>(Z)-ip(l)]-Z + l, (11) 



where ip(n) is the digamma function and £ = f3q 2 /Db is 
the Manning parameter (|] . 



3. Results and Conclusions 



Adding the above contributions we obtain the full 
Helmholtz free energy for the polyelectrolyte plus salt 
system, / = f ent + P c + f cc + /pp. Minimization of the 
free energy allows us to find the number of condensed 
monovalent and divalent counterions, tub and ds- In 
Fig.l we show the dependence of the binding fractions on 
the density of divalent salt, at fixed polyion and mono- 
valent salt concentrations. At low densities of divalent 
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FIG. 2. The osmotic pressure coefficient for the polyelec- 
trolyte multivalent salt system. The parameters have the 
same value as in Fig.l. 



This result could have been anticipated a priori. The 
presence of a condensed divalent counterion does not 
stimulate any further condensation of additional parti- 
cles, but on the contrary, the electrostatic repulsion be- 
tween the like-charged counterions tends to inhibit any 
further association. The situation is quite different in 
the case of binding by the ionic surfactants. After the 
first amphiphilc is associated, the condensation of addi- 
tional molecules is energetically favored since the buildup 
of the hydrocarbon density in the vicinity of a polyion 
helps to exclude water and, thus, reduces the unfavor- 
able hydrophobic energy of the alkyl tails. This, then, 
explains why the association with ionic surfactant and 
lipids exhibits a large degree of cooperativity while the 
binding of divalent counterions proceeds in a completely 
uncooperative fashion. 

Finally, using the number of associated ions, the os- 
motic pressure, defined as (3p = (Zp p + Ipu + ^Pd)4>p, 
can be computed as a function of divalent salt density 
(Fig. 2). 
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